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PULSE -RESPONSE CURVES OF CONVENTIONAL LOW-PASS FILTERS 

By William D. Stanley* and Doris M. Higgins 
Langley Research Center 

SUMMARY 

A catalog of 100 different pulse responses of conventional low-pass filters is pre- 
sented. The pulse responses were obtained by a computer simulation of a general digital 
filter whose characteristics could be varied to correspond very closely with the various 
analog filter functions. 

The standard Butter worth, the Chebyshev (with various ripples), and the Bessel 
filter functions are considered. Filter orders from 1 to 5 are included. For each filter 
type and order, responses corresponding to four different products of bandwidth and 
pulse width are considered. 


INTRODUCTION 

The response of a low -pass filter to a square-pulse excitation is a topic of practical 
interest in the design of systems employing pulse-type waveforms. Practical pulse wave- 
forms employed in radar systems and in coded communications systems may often be 
approximated by the theoretical square-pulse function in many applications. Whenever 
an assumed ideal square pulse is employed in a system, the effects of a given filter on 
the waveform must very often be considered since filtering is required in order to sepa- 
rate the desired signal from undesirable signals and noise. For any finite bandwidth, the 
output of the filter will always be a distorted version of the input pulse. The severity 
and nature of the distortion depend on two factors: (1) the nature of the filter character- 

istic and (2) the bandwidth (or cutoff frequency) of the filter. 

The proper selection of filter types and cutoff frequencies is frequently required in 
order to optimize the design objectives in systems employing pulse waveforms. Usually 
this process necessitates the computation of the pulse responses of several different pro- 
posed designs in order to select the form most suited for the particular application. 
Although the computation of the pulse response is, in theory, a straightforward application 
of fundamental linear-system techniques, the actual details of the calculations may be 
extremely unwieldy for reasonably complex filters. 

^Associate Professor of Engineering, Old Dominion University, Norfolk, Va. 
(Consultant at NASA Langley Research Center.) 
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In this report, a catalog of pulse responses for several different types of low-pass 
filters is presented. For each particular filter characteristic, curves for different 
products of pulse width and bandwidth are presented. The curves were obtained by 
digital -computer computation utilizing a digital -filter simulation. The curves obtained 
should be useful in the design of systems employing waveforms with approximately 
square pulses in conjunction with band-limiting filters. 

GENERAL DISCUSSION 

A comprehensive treatment of the pulse response of the ideal low -pass filter as a 
function of the cutoff frequency and pulse width was provided by Goldman (ref. 1). Fur- 
ther results, which included the effects of various distortion functions, were provided by 
Murakami and Corrington (ref. 2). A more recent treatment of some of the previous 
results was given by Schwartz (ref. 3). 

The results of the mentioned investigations and of similar investigations have been 
employed extensively in analysis of communication systems. However, the major limi- 
tation in these studies is the fact that ideal filter characteristics and ideal mathematical 
distortion functions were assumed. Since the ideal filter is not physically realizable, 
these results are only approximate. In many applications, the error introduced by 
assuming an ideal filter may be very critical. 

Most filter designs today are patterned after the so-called "modern" filter synthe- 
sis approach, in which the transfer function and s -plane concepts are utilized. Notable 
among such filters are (a) the Butterworth, or maximally flat amplitude response, 

(b) the Chebyshev, or equiripple amplitude response, and (c) the Bessel linear phase 
response. Henderson and Kautz (ref. 4) presented a catalog of the step and impulse 
responses of practical low -pass and high -pass filters belonging to the preceding classi- 
fications. Although these curves have been used extensively for calculations of step and 
impulse responses, they are not ideally suited to the determination of pulse responses. 

The present study was made because of a need for a set of pulse responses for the 
conventional filters currently receiving widespread usage. Although several different 
computational approaches could be employed in obtaining these curves, the method actu- 
ally chosen was that of programing a digital -filter simulation of the corresponding analog 
filter characteristic. The bilinear transformation was employed in the approximation 
procedure as described in the appendix. 

The results of the digital -computer simulation were plotted directly on a computer 
software system. As a means of checking for gross error, all the results were quantita- 
tively verified in the laboratory. This task was accomplished by synthesizing the filter 
types with variable components and observing the pulse responses on an oscilloscope. 
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A systematic approach was employed in which the filter characteristics could be readily 
changed by changing the dial settings. A standard pulse generator employing a low-duty 
cycle was used to excite the filters. Although not a precise test, this procedure was 
expected to serve as an extremely useful complement to the computer simulation and to 
assist in uncovering possible mistakes in the basic program. No measurable discrep- 
ancles were observed between the experimental responses and the computer simulation 
results. 


DISCUSSION OF DATA 

The results of this study are presented as 100 different pulse responses in fig- 
ures 1 to 100, A few comments regarding assumptions and curve parameters are in 
order. The input square pulse in all cases is normalized to have unit amplitude and unit 
width. The horizontal scale for the response curves simply reads TIME, The width of 
the normalized input pulse corresponds to 1 on the TIME scale. 

Each filter response curve is characterized by a certain bandwidth B (or equiva- 
lent cutoff frequency). The bandwidth is defined for different filter types as follows: 

(a) For Butterworth and Bessel filters, the bandwidth B is defined as the fre- 
quency at which the response is 3 dB below the de level. Both Butterworth and Bessel 
filters have monotonically decreasing amplitude responses, and there is only one fre- 
quency possessing the 3-dB constraint, 

(b) For Chebyshev filters, the bandwidth B is defined as the highest frequency at 
which the amplitude response is bounded by the prescribed passband ripple. This defini- 
tion permitted extensive usage of the transfer -function -coefficient tables provided by 
Weinberg (ref. 5). 

(c) The first-order low -pass filter can be thought of as the limiting case for all the 
filter types. In this case, the bandwidth B for this filter is the 3-dB frequency in all 
cases. 

The relative bandwidth s for the filter response curves are given as the dimension- 
less product BT. (The bandwidth B is in hertz, or cycles per second, and the pulse 
width T is in seconds.) Since the pulse width is held constant at a normalized unit 
value, the various BT products can be interpreted as different values of the foandwidfchs. 

As an aid in determining the contents and figure number of different response 
curves, a summary of the various filter types and characteristics is provided in table I, 
The decibel specification of a Chebyshev filter refers to the prescribed passband ripple 
for that case. The order of a given filter is the number of low-pass, s-plane poles in 
the transfer function. 
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CONCLUDING REMARKS 

The different pulse responses of conventional low-pass filters have been presented 
in the form of response curves which were obtained by a computer simulation. The filter 
functions considered are the standard Batterworth, the Bessel, and the Chebyshev for 
filter orders from 1 to 5 and for lour different products of bandwidth and pulse width, 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va,, August 29, 1969, 
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APPENDIX A 

DIGITAL -FILTER SIMULATION 
Bilinear Transformation 

The pulse responses presented. in this report were obtained from a digital -filter 
simulation employing the bilinear transformation (ref. 6). The bilinear transformation 
has received a considerable amount of attention recently in the design of digital filters. 
Since all the filter types considered are low pass in nature, with all zeros of transmis- 
sion located at infinity, the Laplace transform prototype p-plane transfer function of 
order k may be represented in the form 

G(p) ^ j7 (AD 

a Q + a^p + a 2 p* + ... + a k p K 
The bilinear transformation is given by 


p= C 


1 - 

1 + z~l 


(A2) 


where C is a mapping constant, and z is the z-transform variable of sampled -data 
theory. 

Substitution of equation (A2) into equation (Al) yields the resulting z -plane transfer 
function G(z). In general, G(z) may be expressed in the form 


b A + b.z-* + b 0 z-2 + ... + b, z~ k 

G(z)= 0 1 2 ~- - k — (A3) 

Cq + CjZ-i + + ... + C k Z“ K 

The various constants in equation (A3) may be related to the constants in equations (Al) 
and (A2) by expansion. A summary of these relationships is given in the following 
section. 

Let x(n) represent the sampled input (pulse function) to the digital filter, and let 
y(n) represent the sampled output (pulse response). The digital computer realization 
of equation (A3) is then given by the difference equation 


y(n) - ^-jb 0 x(n) + b.jx(n - 1) + bgx(n - 2) + ... + b^x^n - kjj 


c 0 



y(n - 1) + c 2 y(n - 2) + ... + c k y(n ~ k) 


(A4) 
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APPENDIX A - Continued 

The mapping constant C is chosen by specifying the correspondence between the 
p -plane and the s -plane. A fundamental relationship of z-transform theory is 

z _ e st s (A5) 

where s is the Laplace transform variable representing the actual system response 
and t s is the time between samples. Let s = jco represent the actual imaginary axis, 
and let p = jX represent the prototype imaginary axis. It can then be verified by sub- 
stituting equation (A5) into equation (A2) that 

Wto 

X = C tan (A6) 

z 

Let X = X r be the prototype reference angular frequency in radians/sec, and let w = u> r 
be the actual reference angular frequency. Assume that the width of the input pulse 
is T, and assume that there are N sample points employed in the interval 0 < t £ T. 
Then, 


T = Nt s (A7) 

Assuming that the number of sample points is chosen sufficiently large to yield a good 
digital -to -analog filter correspondence, the constant C may be determined to a close 
approximation to be 


C = 


NX r 

rf r T 


(A8) 


In all cases considered, the quantity f r was interpreted as the cutoff frequency, or the 
bandwidth B, of the pertinent filter characteristic as explained earlier. With this 
change in notation 


NX r 

?rBT 


(A9) 


where BT may be interpreted as the product of bandwidth and pulse width. 

The computer program was written so that the number of sample points N could 
be varied. It was found that excellent results were obtained in all cases with N = 200, 
This choice results in a Nyquist, or folding, frequency equal to 100/T. Of course, it is 
quite possible that other values of N would have worked as well since no attempt was 
made to determine an optimum value. 
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APPENDIX A - Continued 


Digital Filter Coefficients 

The digital filter coefficients in the transfer functions G(p) and G(z) are as 
follows: 

For k = 1, 

b 0 = a 0 

b l = a 0 

c Q = a Q + a t C 

C 1 = a 0 ~ a l C 

For k = 2, 

b 0 = a 0 
b l = 2a 0 
b 2 = a 0 

c 0 = a Q + a lC + a 2 c2 

= 2Sq - 2agc2 

c 2 = a 0 “ a l C + a 2 c2 

For k = 3, 


b 

b 

b 

b 

c 


0 = a 0 
1 = 3a 0 
2 = 3a 0 


3 - a Q 

0 = a Q a 2 C + agC 
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APPENDIX A - Continued 


Cj = 3a 0 + ajC - a g C 2 - 3a 3 C 3 
c 2 = 3a 0 - a^C - agC 2 + 3a,jC 3 
c 3 = a Q - ajC + agC 2 - a 3 C 3 

For k = 4, 


For k = 5, 



c q a Q + ajC + a 3 C 2 + a^C 3 + a^C 4 
cj = 4a 0 4- 2 ai C - 2a 3 C 3 - 4a 4 C 4 

c 2 = ® a 0 “ 2a 2 < -'^ + ^a^C 4 

Cg - 4a Q - 2ajC + 2agC 3 - 4a 4 C 4 

^4 — ag - a^C + agC 2 - a 3 C 3 + a 4 C 4 
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APPENDIX A - Continued 


b 5 = a 0 

c 0 = a.Q + a^C + a 2 C 2 + a^C 3 + a^C 4 + a g C 3 
= 5a Q + 3ajC + a 2 C 2 - a g C 3 - 3a^C 4 - 5a g C 5 
c 2 = 10a Q + 2a t C - 2a 2 C 2 - 2agC 3 + 2a 4 C 4 + 10a g C 5 
eg = lOap - 2a^C - 2a 2 C 2 + 2a g C 3 + 2a 4 C 4 - 10a g C^ 
c 4 ■■ 5arj - 3ajC + a 2 C 2 agC 3 ~ 3a 4 C 4 + 5a g C^ 
c 5 = a Q - a^C + a 2 C 2 - agC 3 + a^C 4 - a g C^ 


Analog Filter Coefficients 

Detailed consideration of filter characteristics are provided in such synthesis 
texts as Weinberg (ref. 5) and Van Valkenburg (ref. 7). The analog filter coefficients 
are given in the following table: 


Filter 

a 0 

a l 

a 2 

a 3 


k = 

1 



All 

1 

1 




k : 

2 



Buiterw or til 

1.0000000 

1.4142136 

1.0000000 


Chebyshev (0.5 dB) 

1.5162026 

1.4256245 

1.0000000 


Chebyshev (1 dB) 

1.1025103 

1.0977343 

1.0000000 ; 


Chebyshev (2 dB) 

.6367681 

.8038164 

1.0000000 


Chebyshev (3 dB) 

.7079478 

.6448996 

1.0000000 


Bessel 

3.0000000 

3.0000000 

1.0000000 



k = 

= 3 



Butterworth 

1.0000000 

2.0000000 

2.0000000 

1.0000000 

Chebyshev (0.5 dB) 

.7156938 

1.5348954 

1.2529130 

1.0000000 

Chebyshev (1 dB) 

.4913067 

1.2384092 

.9883412 

1.0000000 

Chebyshev (2 dB) 

.3268901 

1.0221903 

.7378216 

1.0000000 

Chebyshev (3 dB) 

.2505943 

.9283480 

.5972404 

1.0000000 

Bessel 

15.0000000 

15.0000000 

6.0000000 

1.0000000 
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APPENDIX A ~ Concluded 


Filter 

a 0 

a l 

k 

El 

a 3 

a 4 

I a 5 I 

Butterworth 

1.0000000 

2.6131259 

3.4142136 

2.6131259 

1.0000000 


Chebyshev {0.5 dB) 

.3790506 

1.0254553 

1.7168662 

1.1973856 

1.0000000 


Chebyshev (1 dB) 

.2756276 

.7426194 

1.4539248 

.9528114 

1.0000000 


Chebyshev (2 dB) 

,2057651 

.5167981 

1.2564819 

.7162150 

1.0000000 


Chebyshev (3 dB) 

.1769869 

.4047679 

1.1691176 

.5815799 

1.0000000 


Bessel 

105.0000000 

; 

105.0000000 
1 r 

45.0000000 
- ^ 

10.0000000 

1.0000000 


Butterworth 

1.0000000 

K 

3.2360680 

— □ 

5.2360680 

5.2360680 

3.2360680 

1.0000000 

Chebyshev (0,5 dB) 

.1789234 

.7525181 

1.3095747 

1.9373675 

1.1724909 

1.0000000 

Chebyshev (1 dB) 

.1228267 

.5805342 

.9743961 

1.6888160 

.9368201 

1.0000000 

Chebyshev {2 dB) 

.0817225 

.4593491 

.6934770 

1.4995433 

.7064606 

1.0000000 

Chebyshev (3 dB) 

.0626391 

.4079421 

.5488626 

1.4149847 

.5744296 

1.0000000 

Bessel 

945.0000000 

945.0000000 

420.0000000 

105.0000000 

15.0000000 

1.0000000 
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TABLE I.- PRESENTATION OF RESULTS 


Filter type 


First-order 


Butterworth 


Chebyshev (0.5 dB) 


Chebyshev (1 dB) 


Chebyshev (2 dB) 


Chebyshev (3 dB) 


Bessel 


Butterworth 


Chebyshev (0.5 dB) 


Chebyshev (1 dB) 


Order 


BT, Hz-sec 


Figure 
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TABLE I.- PRESENTATION OF RESULTS - Continued 


Filter type 
Chebyshev (2 dB) 

Chebyshev (3 dB) 

Bessel 

Butterworth 

Chebyshev {0.5 dB) 

Chebyshev (1 dB) 

Chebyshev {2 dB) 

Chebyshev (3 dB) 

Bessel 

Butterworth 
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Figure 3.™ Response of first-order filter tv it h BT “ 2, 
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Figure IL- Response of second-order Chebyshev (0.5-dB rippled filter with BT = 2. 
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Figure 19.- Response of second-order Chebyshev (2-dB ripple! filter with BT - 2. 












-# illl l l l U 


JJ I f II 1 1 1 1 


l.Q- 


JJJJJJ i f I 


i.B 


_n 1 1 n r n h 1 1 1 m n 




time 


MNIllii 


3,0 


hi i ! 1 1 f r 


3*5 


1 1 ii m i j I 


* 1*0 


i ii i r m 1 1 


4,5 


mnimJ 


5 .Q 


Figure 26^ Response of second-order Bessel filter with BT - t 
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Figure 30,“ Response of third-order Botterworth f liter with 81 “ L 
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Figure 31- Response of third -order Butterworfh fitter with 81 = 2 , 
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Figure 32 Response of third-order Butterworfh filter with BT = 5 
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Figure 40.- Response of third-order Chebyshev ( i-d B rippfef filter with BT = 5. 
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Figure 43.- Response of third-order Chebyshev I2-6& ripple) filter with 81 - 2 t 
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Figure 53.- Response at fourth -order Butterworth filter with 8T - 0.5. 
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Figure 76,- Response of fourth-order Bessel filter with BT = 5, 
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Figure 84.^ Response of fiffh onder Chebyshev (0.5-dB ripple) filter with BT = 5, 
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th-ortfer Chebyshev (2-dB rippiet fitter with 8T = 5. 
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